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Introduction 



In this note, we consider the following process^ with stochastic volatility : 

dF{t) = a{tYc{F{t)) dW{t) 

dcr(t) = X'in' - a{t))dt + va[tY dZ{t) . (1) 

¥.{AW{t)dZ{t)) = pAt 

• for 6^ I, c{F) ^ and A' = we get the SABR model; 

• for (5=1/2 and c{F) — F we get Heston's model. 

A direct application of Feller's criterion^, detailed in appendix 1, shows that for 1/2 < 5 ^ 1 there 
is no explosion for the stochastic volatility, that is to say 

p[{yt ^ 0) {(T{t) e]o,+oo[)] = 1. 

This result remains true for S = 1/2 if we add the condition 2A/i/i'^ ^ 1. We suppose in the following 
that we are in such cases. 

Let us write Gf.t.{t, f ,a) he the density of probability to get to the point {F,Yi), leaving from the 
point {f,cr) after a time r. Then G follows the classical differential equation (where A' and fi' .should 
be adjusted!^ ) 

The goal for this note is to get suitable approximations for the 6-model, for small t. There are three 
steps to get them. 

• We first explain how to convert our stochastic problem into a geometrical one. We introduce 
the Laplace Beltrami operator to have a geometric point of view, and then use a suitable isometry 
to work in a simpler space, called the S -space. 

• Thanks to this geometric point of view, we then give a first-order asymptotics for the prob- 
ability density G. Moreover, we want this expression to be computed for all F and d without 
having to calculate any integral. The main tool for this part is Molchanov 's theorem^ . 

• The third one is to get a first order (in t) approximation for : 

^In the first chapter, with generahties about stochastic volatility problems, wc will understand why we consider this 
particular process 

■^see lOANNlS Karatzas, Steven E. Shreve, Brownian Motion and Stochastic Calculus, Springer. 
^ Jim Gatheral, Stochastic Volatility and Local Volatility, Case Studies in Financial Modelling Course Notes, Courant 
Institute of Mathematical Sciences, Fall Term, 2004. 

^S. A. Molchanov, Diffusion Processes and Riemannian Geometry, Russian Math. Surveys 30 : 1 (1975). 
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— the transition probability from (cr, /) to F, where we integrate over E because the final 
volatility is not observed : 



This gives us the price for digital options. 

~ the implied volatility (735 (A'), defined by the equality of the prices of calls for (1) and a 
Black-Scholes model with volatility obs{K) : 



Such approximations are given in literature for the SABR model, that is to say 5 = \. The most 
successful closed form formula giving an approximation for options is Hagan's"^ . In a preprint^ Andrew 
Lesniewski uses a geometrical method (in the hyperbolic plane) to get the density's asymptotics, but he 
uses a Dyson formula instead of Molchanov's theorem. Consequently our results are sensibly different. 

Berestyckv' gives an evolution equation for the smile. His PDF is generally difficult to solve and 
hides the geometrical aspect of the problem. This is the reason why, for our 5 -model, we do not try to 
use the available PDF and directly calculate the first order. 



^Patrick S. Hagan and Diana E. Woodward, Equivalent Black Volatilities, Applied Mathematical Finance 6, 
147-157 (1999). 

'^Patrick S. Hagan, Deep Kumar, Andrew S. Lesniewski and Diana E. Woodward, Managing Smile Risk. 
^Patrick Hagan, Andrew Lesniewski, Diana Woodward, Probability Distribution in the SABR Model of Stochas- 
tic Volatility. 

^Henri Berestycki, Jerome Busca and Igor Florent, Computing the Implied Volatility, Communications on 
Pure and Applied Mathematics, Vol. LVII, 0001-0022 (2004) ©2004 Wiley Periodicals, Inc. 




,tie;(cr, K). 




^dT 



A geometrical view of stochastic volatility 

problems 



We consider the following general process, 

d^^ = /i(E,F)dt + f2{^,F)dWi 
dE ^ gi(S)di + g2(S)dVF2 ' 

with < dW^i, dW2 >= pAt. We call G the transition probability density, from a state (/o, <Jq) to a state 
(/, a) . Then the backward Fokker Planck equation can be written 

9rG - i {fidffG + gldaaG + 2phg2d.fG) + /iS/G + gid^G. 
The Laplace Beltrami operator is defined as 



for a metric G. If we consider the upper half-plane with metric 



G 



1 



then we have 



gl 



-pf2g2 
/I 



^1 
p/252 



/I p/2g2 



.g = detG = 



1 



where / is a first-order operator. In order to apply Molchanov's theorem^, we need to be able to 
calculate the geodesies of the space with metric 



G = 



(1 



^ h 92 
, 2 

1 



92 



The purpose of this chapter is to understand the sufficient conditions for the geodesies to be easy 
to calculate. 



1.1 The isometry 



In order to be able to calculate the geodesies easily, we look for a suitable isometry. We impose the 
E-coordinate to be stable (simpler case to keep the upper half-plane stable) , so we look for an isometry 
on the form 

1 -p 







1 



^and more generally Minaksliisundaram-Pleijel recursion formula. 
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CHAPTER 1. A GEOMETRICAL VIEW OF STOCHASTIC VOLATILITY PROBLEMS 



Then H = (V$)-i G(V$)-i is diagonal iff 



92U = /2W. 



We then have 



H = 



J 2 








92 



For <J> to be well defined, the following condition is necessary : 

d^u = dpv. 

Interesting special cases are for u = f2 and, more generally for M//2 not dependant on F. 

1.2 Special cases 

In this part, we consider special cases showing the interest of this geometrical formulation. The first 
one confirms results about non-stochastic volatility models. The second one gives a condition on /2 
and 92 for both horizontal and vertical lines to be geodesies, which is a very useful case. Finally, the 
third one shows that the (5-model is the only « /2-separable » model for which we have the homogeneity 
property. 

• .92 = (no stochastic volatility) : then the geodesical distance between {fo,ao) and (/, cr) is 
infinite, except if ctq = cr, so the probability of transition is concentrated on this line of the 
plane : cr = cto- Varadhan's theorem in dimension one then gives 



log(P(/o ^ /)) ^2:, with d 



dx 



Ifo f2{cro,x) 



we want to find an isometry between G and a simpler space H. The most interesting is a function 
$ such as 



V$ = 



VT^32 

1 



because we then have 



G = (V<I>)^iJ(V$) with H 



However, to find a convenient function $, the partial derivatives must follow Schwarz's theorem 
of commutation (Poincarc's criterium for the differential to be exact). This is the case iff /2 does 
not depend on S. So this case is simply not pertinent^. 

if u/ f2 only depends on E, then with Schwarz's relation one can easily check that the function 
/2 is separable : let us write /2 = a{F)b(T,). Then we also want to find an isometry between G 
and a simpler space H. The good function $ is such as 



V$ = 



1 



^However, one can imagine new models where 32 depends on F, which seems quite natural. To be able to calculate the 
geodesies, we then need the relation 9e(1//2) = p/52) (*), which just necessitates well chosen functions. Moreover, 

the relation (*) gives interesting qualitative results : if the correlation is positive, then the volvol decreases with F ; the 
contrary if the correlation is negative. 



1.3. SUMMARY 
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because we then have 



G = (V4>)^i?(V$) with H 









1 

if 



Schwarz's commutation relation is now true, so we really defined an isometry. We now just need 
to be able to calculate geodesies in the space H. First notice that the metric does not depend on 
the coordinate x, so, using Killing Fields, a system of equations for the geodesies is 



y 



Cl 
C2 



with {x,y) = E) 



Thanks to this system, one can check that wc have the homogeneity 



property (if 7 is a geodesic, so is A7) iff b{y) = .92(2/) = on the upper half-plane"^ 



1.3 Summary 



Proposition. We consider all isometries of the upper half-plane which are identity for the a-coordinate. 
Then : 

• if we want the resulting metric to be diagonal and only T,-dependant, then /2 must be separable 
(f2 — a{F)b{Y,)). One then can effectively find a suitable isometry and the resulting metric is 
1 

// the geodesies follow the homogeneity property, then the only 



given by H - 
suitable metric is H 



6(S)2 





1 






1 

^27 



one can find other isometries giving the following simpler metric : H 



1 
4 



This ne- 



cessitates the dependance 32(/jf) (so this is not exactly a stochastic volatility model) and the 
equation d^{l/f2) = dF{-plg2)- 



•^Indeed, we have such homogeneity for every geodesic iff the following system is satisfied for every A and y in (we 
also suppose that the functions b and 32 are strictly positive and '^^^ on K^) : 

b{\y) ^ b{y) 

b{\y) = c{\)b{y) 

This implies that we also have g2{Xy) = c{X)g2(y), so we only need to show that 6 is a power. Moreover, as b{\y) = 
c{\)b(y) = c(j/)6(A), we have b{y)/c{y) = est, so we only need to show that c is a power. 

For n e N*, we have c(A)"6(?/) = c(A"-l)6(A?/) = ■ ■ • = c(A)6(A"-ly) = b(A"y) = c{\")b(y), so c(A") = c(A)". The 
following is straightforward : we also have c(A'") = c(A)'" for r £ Q*^, so, as is dense in R_|_ and c is continuous, we 
have c(A') = c(A)= for s eR+. U y > 1, we then have c(y) = c(e'"!') = c(e)'"!' = y^'^^'-'^l = y^^ . For < y < 1, we also 
have c(y) = c{(l/e)-l"f) = c{l/e)-l"f = y-ln^^Cl/e) = yS2 , As c is at y = 1, we have Si = fo, so c, b and 92 can be 
written y^ . 



CHAPTER 1. A GEOMETRICAL VIEW OF STOCHASTIC VOLATILITY PROBLEMS 



The (5-model : geometrical formulation 



In order to have a geometrical vision of the situation, we first change the variables and the function 
(approximately like in Lesniewski's preprint^), writing : 

• s = v^T ; 

• a; = /; 

• X = F; 

Now the Green function, solution of the PDE (2), is solution of the following system (where we 
have A = W jv^ and /i = ^ jv) : 



OK 
ds 



1,. 25-2, ,2.5 

2'^ y 



2pc{x) 



dxdy dxdy 



Kx,Y{Q,x,y) 



5{x,X) 5{y,Y) 



+ A(y-M)f 



(1.1) 



In this part, we first rewrite the PDE (1.1) in terms of the Laplace Beltrami operator of a suit- 
able space (the important properties of this operator are summarized in appendix 2). We then use an 
isometry that allows us to consider a simpler generic space. Our results are summarized in a final table. 



2.1 Suitable form for the partial differential equation 



The metric associated with the PDE (1.1) is 

1 



G 



1 —pc{x) 
-pc{x) c{x)'^ 



We call the corresponding Riemannian space G. We also need the determinant 5 of G and the 
inverse matrix : 



1 



9 = 



The Laplace Beltrami operator is defined as 



c{x)'^ pc{x) 
pc{x) 1 



-^Andrew Lesniewki, Probability Distribution in the SABR Model of Stochastic Volatility, preprint 
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CHAPTER 2. THE S-MODEL : GEOMETRICAL FORMULATION 



A simple calculation gives A'^ = v'^^-'^y^\c{xfd^x + 2pc{x)dxy + dyy) + iy'^^~'^y'^^c{x)c'{x)dx. This is 
the reason why (1.1), with the suitable normalized initial condition, can be written as 



with 



Xf(0,z) = S{z,Z) 



(1.2) 



(1.3) 



and (5(z, Z) = 5{x, X) S{y, Y). The interest of this formulation is that wc have separated the « intrin- 
sic » Laplacian from the first-order term (without changing variables like in Hagan's papcr^). 

2.2 The isometry 

Wc follow here the general methodology of the chapter 1. Let us introduce H, the (5-spacc with metric 

1 



H 



-4 



Let us $ be the application of the upper half-plane 

1 



X 

y 



V 



where p is a positive constant. Then the jacobian matrix of <i> is 

1 p 



1 



,8-1 



1 



so one can easily check that between the metrics the following relationship holds : 

G = (V$)^iJV<I>. 

This relationship exactly means that $ conserves the scalar product : in G, for two trajectories 

. T 

7^ and 7^ crossing in z at t, < d/dt <I>(7^),d/di $(7^) >h= {V^{z)-i{Y H{V^{z)-il) = 7^ G-fl =< 
d/d< 7f,d/di 73 >g. This local property implies the global result that <f> induces an isometry between 
G and H. 

The following paragraphs are quite technical : in order to simplify our problem, we find a diffusion 
equation in H equivalent to the one in G. We then find the relations between the transition probabilities 
associated to these diffusion equations. 



2.3 Reformulation of the hyperboHc problem 



Wc now consider the following Green problem in 



KtJ^{^,z') 



' 5{z'MZ)) 



det [V^{Z)]~^ 5{z,Z) 



(1.4) 



with 



^Patrick S. Hagan, Deep Kumar, Andrew S. Lesniewski and Diana E. Woodward, Managing Smile Risk. 



2.4. FROM DENSITY TO PROBABILITY OF TRANSITION 13 



• z' = $(z) ; 

• the Laplace Beltrami operator corresponding to the Riemannian space H, therefore A''^ = 

• = V$($-i(v4))/^($-i(v4)). 

If K"^ is a solution of (1.4), then det [V$(Z)] o <i> is a solution of (1.3) (det [V*(Z)], which is a 
normalization constant, is written C in the following). Indeed, as <f> does not modify the time s, 

dK^ _ d{C A'H o $) 
ds ds 

- C Qa^^h^ o$ + C V^/^IX'^) 
= C ^A'^{K^o<P) + C f{K'^o<S>) 
dK^ = QA« + /«)if«. 

Thanks to the normalization constant det [V$(Z)], follows the condition at s ~ 0, so from the 
uniqueness of the solution of (1.3) we deduce that 

/C^ det [V$(Z)] A'^ o (1.5) 

2.4 From density to probability of transition 

The quantity K'^ is a density with respect to the euclidean metric. To be consistant with the next 
chapter, we now note K the corresponding density probability with respect to the metric H : 

K = = y^'K"^. 



We have then 



Vdetff 
GdS 

K^{X^,Y^)dY^ 

o $(X*^, Y^) detV$(X'^, Y^) dY^ 

^° '^ff'^''^ detva>(x«,y«) dy"^ 
A-o$(x«,y«) 1 

yH2<5 ^2i-2yj-r^c(XG) 
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CHAPTER 2. THE S-MODEL : GEOMETRICAL FORMULATION 



where D is the image of the straight hne x = X"^ by the appUcation <I>, that is to say the straight hne 
with equation 



du 



So if wc define 



then 



P 



1 



(1.6) 



In the following wc will need the expression of the moment of order 25 of G, that is to say 



/r 



If we define 



then the same calculation as above gives 



M 



Here is a table that summarizes the results of the first part 



(1.7) 
(1.8) 



From G to 



The space 
The metric 



1 —pc{x) 
-pc{x) c{x)'' 



H -■ 



The isometry 



The LBO 







(1:1 





The first-order operator 



1 2(5-2 25 ^ \ // N 
-it/ y c{x)c (x) 

My - m) 



f*{A) =V'S>{'l>-\A))f'{'S.-\A)) 



The diffusion equations 
Link between K'^'' and 



5{z,Z) 



' 5(^',<I>(Z)) 



X*'' = det [V$(Z)] if" o * 



The transition probability 



P = 



pH _ f K{X^, 
JB y,H2 



The moment 



M 



All relationships above imply that we just need to get a suitable approximation for K''^. This is the 
purpose of the next chapter. 



From Molchanov's theorem to the 5-model 



We have shown in the previous part that it is sufficient to calculate all probability densities in H, the 
(5-space. This can be done thanks to the following theorem, expressed in terms of geometric invariants. 
The necessary calculations and adaptations to the (5— space will be done after. 



3.1 Molchanov's theorem : the 0-order expansion 

Molchanov's theorem gives an asymptotic behavior for the density transition function, for small ma- 
turities. More precisely, we use the following notations : 

• I is the dimension of the space ; 

9 H is a specified metric of a space H ; 

• is the Laplace-Beltrami operator corresponding to the metric H ; 

• / is a first-order operator ; 

• K{s, a, b) is the density function (with respect to the clement of volume v^) corresponding to 
the diffusion process ^ = (jA"^ + /) K with initial condition K{Q^a,b) = (5(a, 6) : K{s,a,b) 
represents the probability density of getting in b, leaving from a, after a time s, for the metric 
H. 



Supposing sufficient conditions of regularity (uniqueness of the geodesic between two points. . .), we 
can now formulate the theorem we need. 



Molchanov's theorem. Under all the conditions above, if d{a, b) is the distance between a 
and b along the unique geodesic linking them, j, then 



Kis a b) ^_^-ii^+wia.b) dia,by~^ 



with 



• 4'(a, 5) = dS/ dip where a cone of light from x with a solid angle d^} illuminates an area 
dS on the geodesic hypersurface orthogonal to j at y (cf the following figure ) ; 

• W{a, b) ~ Jq^"^'^"^ < d'y / ds > ds : this is the work of the field f along the geodesic. 



Proof. A complete proof can be found in Molchanov's article-'^. However, one can give an intuition 
of this result if there is no first-order operator, by analogy with the euclidean case. 



^S. A. MOLCHANOV, Diffusion Processes and Riemannian Geometry, Russian Math. Surveys 30 : 1 (1975). 
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16 CHAPTER 3. FROM MOLCHANOV'S THEOREM TO THE S-MODEL 



For a flat space, the above formula is clearly 
true, and for a curved space the density natu- 
rally decreases with dS/dup, where d(p is the ini- 
tial angle between geodesies leaving O. Indeed, 
imagine a diffusion in this space, materialized 
with particles. The particles are diffused in an 
isotropic way from O, and they remain blocked 
between the geodesies, so KdS is analogous to 
a one-dimension diffusion density with a density 
Kod(p in O. 

This does not leads us to the exact formula 
(because particles come from the other side of 
dip), but it gives an intuition of the origin of the 
function ^. 




Figure 3.1: Geodesic surface for a curved space 



From a computational point of view, we need a definition of ^ relatively to coordinates. Let (e^) 
be a base of the orthogonal to 7 in O. Let us ^i{t) (l<i<Z— l,l<t< d(a, h)) be the Jacobi field^ 
along 7 such as = and /[{O) = a, and ^{t) = [Ji{t), . . . , /i-i{t)]. Then 

^{x,y)= deiSf {d{a,b)). 

Indeed, det [/i, . . . , fi-i] represents the area engendered by a polyhcdra of sides /i, . . . , so 
dei2f{t) represents the infinitesimal area engendered by the geodesies at time t, fo a solid angle 
dp. □ 

What we will need in reality is both of the following corollaries. 



Corollary 1 : dimension 2 (for the ^-model). In the plane, we consider the following 
system, 

Z + R{t)Z = 
ZiO) = , 
Z{0) = 1 

where R(t) is the Gauss curvature for the point j{t) of the geodesic linking a and b. Then 

1 



K{s,a,b) ~ e 



-+^■(0,6) 



' d{a,b) 
Z{d{a,b))' 



Proof. Wc only need to show that the the evolution of the length Z of the unique vector of 3f is 
governed by the equation above. The evolution equation for a Jacobi field is 

y'' + y^{t)R[e,{t),^,^,ek{t)]=0 
for each of the I — 1 composants j/j, , where : 

• i? is the Riemann tensor ; 

• {ei(0), . . . , e;_i(0)} is a base of the orthogonal hyperplanc of 7 at a ; 

• e,;(t) is the parallel transport of 6^(0) along 7 relatively to the field 7 ; 



^see JuRGBN JOST, Riemannian Geometry and Geometric Analysis, Springer Univcrsitcxt, 2002. 



3.2. THE FIRST-ORDER EXPANSION 
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• yi is the composant of the Jacobi field on e^. 

We also know that the parallel transport relatively to 7 transforms ei(0) in the orthogonal of 7, 
with norm I7I, written 7-^-. So the evolution equation can be written, in dimension 2. 

y + yR[l'^, 1,1,1'^] = 0. 

From the definition of the Gauss curvature, i?[7-'-, 7, 7, 7-'-] = R{t), which completes the proof. □ 



Corollary 2 : dimension 2 and constant Gauss curvature (for the (5-model v^fith 
^ = 1 : SABR). // the plane has a constant gaussian curvature —k < 0, then 



K(s,a,b) ^ ^,-^+w(a,.) / Vkd{a,b) 
' ' s~.o 2tts y sinh Vkd{a,b) 



Proof. Immediate from the corollary 1. □ 



3.2 The first-order expansion 

We can get a better approximation for the probability density K. Indeed, one can show that for a 
sufficiently regular diffusion and space, K admits the following expansion'^, for any n e N : 



K{s, a, b) = ^^-L_e-^+^(<^'^) y'^|2^(l + /^(a, b)s + ■ ■ ■ + I^a, b)s- + 0(s"+i)). 

Our purpose is to give the equation verified by Ki in the case of the 5— geometry, and then evaluate 
its solution for a point b near from a. 

The expansion above can be written at first order as 

K ^-e-'&g{l + Kis + 0{s^)). (2.1) 
s 

Injecting (2.1) in the diffusion equation dgK = (l/2 + /) K, we successively get the following 
diffusion equations satisfied by the functions d, g and Ki, for the diffusion in the 5— space : 

y^'lVdp = 1, (2.2) 

g = 2dg(^^A'^ + f^d + y^'Wg-Wd\ (2.3) 

1 . m „\ 2(5 



gK, = {-A'' + f]g + y''dgVK,-Vg. (2.4) 



Here | |, • and V reference respectively the euclidean norm, the euclidean scalar product and the 

g), and A« = y^'{d^,,,+dyy) 

(2.2) just traduces that ds^ = {dx^ + dy'^)/y^^ ; 



euclidean gradient (V = {dx,dy)), and A"* = y^^[dxx + dyy)- These equations need some comments 



^see Elton P. Hsu, Stochastic Analysis on Manifolds, American Mathematical Society vol. 38, Graduate Studies in 
Mathematics, 2001, theorem 5.1.1 p. 130. 
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• (2.3) is a transport equation that can be solved by integration along the geodesic : the explicit 
solution is given by Molchanov's theorem : 



9ia,b) ^ —e 



1 „w(a,fc) / d{a,b) 



27r y Z{d{a,b))' 

• (2.4) is also a transport equation, that can be solved by integration along lines directed by Vg 
instead of Vd : the knowledge of the geodesies is not sufficient anymore. This is the reason why 
we will just calculate its initial value {Ki(a, a) = (/^i)o), which will be a sufficient correction for 
us in the following. 

A great calculation (see appendix 2.3) gives us the following result for (i^i)o : 
t,^^ ^ (div/)o (/y)o^3|/oP 

The main result of this part is then the following expansion, in the (5-space : 



with : 

• the functions d and Z computed in 0(1), thanks to the homogeneity property of the (5-space (cf. 
appendix 2.4) ; 

• the function W easily computed as the integral of a known function along a known path. 

• (A'i)o given by the formula above. 

Remark 1. When / = 0, our result for (i^i)o is a special case of a general formula : in any 
riemannian space we have 

(i^i)o - |, 

where R is the scalar curvature. 

Remark 2. A formula by Minakshisundaram and Pleijel'* gives a general asymptotic expansion 
for K to any order, but it is too complex and non effective for the second order. 



^S. Minakshisundaram and Pleliel, Some properties of the eigenfunctions of the Laplace- operator on Riemannian 
manifolds, Canadian J. Math. 1 (1949), 242-256. 



Applications of the previous asymptotics 



In this chapter, we integrate the asymptotic expression obtained previously. This leads directly to 
digital option prices. 

Then, using Hagan's local volatility formula^, we compute, through equivalent local volatility, 
vanilla option prices. 



4.1 The Laplace transformation 

We want to evaluate and defined by (1.6) and (1.8). We use here a Laplace expansion of order 
2 : for / and g sufficiently regular functions, and d with minimum at uq and non-zero second derivative, 
we have 



e » {f{u)+sg{u))du 



2-KS 



fiuo) 



V $"(^io) ^ 
f(^o)$(3)(^t„) 
2$"(?/o)2 



9{uo) + 



f"{uo) <i>W(^o)/(uo) 



24$"(iio)^ 



8$"(?/o)2 



minimal. This is done in appendix 2.5. 

Unfortunately, we are unable to 
calculate the third and fourth varia- 
tions of $ for any strike, but this can 
be done for K = /o, at the money ; 
including these results in the 1-order 
term above will be sufficient : the ef- 
fect of the 1-ordcr correction is more 
important a the money. 

To sum up, we will take the follow- 
ing approximation for P : 



2t:s 



^"{uo)\ 



K 



/(uo)|^+.5 (Kl)0 




/"(^o)| 



K=h 



$(4)( 



Wo, 



\K=fo 



2^"{uo)\ 



K=fo 



8$"(uo)2| 



/'(uo)|^=-'«$(3)(yQ)|if=/o 5$(3)(yQ)2|A:=/o 



2$"(uo)2|^=/o 



+ 



24$"(wo)3|^=/o 



with the following formulas : 



^Patrick S. Hagan and Diana E. Woodward, Equivalent Black Volatilities, Applied Mathematical Finance 6, 
147-157 (1999). 
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urations 



■^d{Yrmn) = ( ^ j^^^L; ± ^FfT COS 6*1 1 , where the sign depends on the foUowing config- 




This second variation for the 5-model is the most important resuh of this part because it gives 
the leading order for pricing the digital options. 

• = dd" ; 



<j,(4)|/^=/o „ -^(4+7^) 



These results are justified in appendix 2.5. 



4.2 The transition probability : pricing digital options 

All the previous results give a methodology to calculate the transition probability P for any S. For 
instance, in the case 6=1, the explicit calculations are done at the end of this paragraph. What 
is done here for ^ = 1 can be done for any 6 : the trajectories of 5-geometry and the Jacobi fields 
calculations are given in appendix. 



4.2.1 Calculation of the field and its work 

If c{x) = (0 < /? < 1 : the case (3 = \., needed for Heston's model, can be handled apart, the same 
way), the function $ has the following inverse. 

Moreover, 

V<i>(7'^(t))/(7«(i)) 



(i-/3)(vT~^x«i + pyi 

yH 



1 



-iy«''c(x«)c'(x«) 



yG2^^^G/3-l _ 



r(7"(0) 



\^2(l-/3)(0^XH+pyf«) 

A {Y^ - /.) 



pA(y^-M) 



From (5- geometry, we know that the hyperbolic space (ie the case 5=1) has constant Gauss 
curvature : this specificity makes all its interest, we can apply the corollary 2. We know that geodesies 
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are circles, so Y'^iX) = y (X% - X^f^^)^ + Y% - {X^ - Xf^n We can then calculate explicitly 
the work W of on these circles. Simplifying by dt, we get 



j JjU 



7H \^2(l - /3) (v^i"^^^ + 



2v/l~p2(l - /3) yr^^H + pyH v^l"^ A- 



This expression does not depend of the choice of time, which is normal for the work of a field on a 
trajectory. 



4.2.2 The 1-order expression of P : the SABR case 



We first need a primitive (in x) of the function 



a + \/l — p^x + pVl — x'^ 



We will note it as G(a, p, x). An explicit formula for G is 



G{a,p,x) 



1 



4Va2 - 1 



-8ap(3 - + a^(-3 + V)) arctan 



- - 1 ( Vl - yo^ + 4ax - Sap^a: - 2i/l - p'^x^ + 4(1 - 4a^)/ci^ arccosx 



2p((--4ay/l-p^+x) v^l - a;2 + (-3 + Ga^) 



arccos x 



We then have 



1 



0-^3/2^2^ 
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where the first-order term has been calculated thanks to the previous variations. The notations are 
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given in the following. 

C = 



Yo 



I 



a 1-/3 

- 2pC + 1 
1 ff'~^ 



y 1-/3 



1 /f^-" 



(^.„d + y.„d)-(^o +yy) 

2(X,„d - Xo) 



/(<?) = ^{iVi- 



q-' + arcsin g J 



wo = r G a, p, — G a, «, 

2(l-/3)^l-p2 V V »■ 7 V 

WRevert = —n ^ r f — f 

v/l-p2 V V J 



d{a^b) — lo, 



i + c - p 



1 - p 



4.2.3 Graphics 

In the SABR case, here are the evolution of the transition probability and an example of a digital 
option that we got. 




The 1-ordcr expression for P seems to be relevant approximatively for r G [0,3], as shown in the 
following graphic : for t > 3, the 0-order expression is more stable : it remains aproximatively a 
distribution function, whereas the first-order does not. 
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Figure 4.1: 1 — J^P for the 0-order (up) and the 1-order (down) of P 

4.3 The implied volatility 

4.3.1 The local volatility : first order asymptotics 

Local volatility can be defined as 

aK{T, f, afdT = E ({dF{T)f \ F{T) = K 
In the case of the (5-model, (1) implies that 

aK{T, /, af = C{Kf¥. ({a{T)f \ F{T) = K 



If we introduce the transition function Gp,-s{t, f^cr), then a simple consideration such as P{A \ 
B) = P{A n B)/ P{B) (except that we reason here with densities and zero-probability events) shows 
that 

Jo Gk,e(t", /,cr)dE Pk{tJ,(t) 

where 

• Pk{t, /, a) = Gk,t.{t, f, a)dTj is the marginal probabiHty distribution ; 

• AIk{t, f, a) = 'E,^Gk,t,{t, /, cr)dS] is the conditional moment of order 26. 

At the leading order Pk and Af^^ have almost the same asymptotic expansion (as we can see with 
a saddle point method) : 

Ml,{Tj,a) ^ j:^i^PK{Tj,a), (2.2) 

T— ►0 

with : 

m y - - uY'^ ■ 

• Y'^ = ,,<5-lyH . 

mm mm ' 
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• y^j„ the point of the straight hne with equation 

v^-^y/Y~^ \Jp c{u) y) 

minimizing the distance to the point $(/, c/i^). 
We can conclude with (1) and (2) that 

aKir,f,a)^i,.'YZJciK)+oir). 

Unfortunately, we don't have a closed form 
for the distance from a point to a straight line 
in any i5-space. However, we can express in 
terms of a new 2-arguments function. 

For this, imagine the following problem. In a 
5-spacc, a rocket is thrown perpendicularly to a 
straight line and crashes on another one. Then 
we define 

16(01,02) = -. 

The function g is well-defined because if 7 is a 
geodesic, so is A7 

If we take 61 and 62 such as 

• tan6'i = — -y/l — p^/p, 

• ta.iid2 = ^/T~f^/{v/a du/c(M) - p). 

Figure 4.2: The function is{0i,02) = 2/1/2/2 

then 

(Jk{tJ,o) = {ais{6i,62)f c{K) + 0{t). 
li V = Q (no stochastic volatility), one can check that we have the expected value 

<7k[t, f,a)=a'c{K) + OiT). 




4.3.2 From local to implied volatility 

The previous calculation made it possible to transform the original stochastic volatility problem in a 
local volatility problem. This is the reason why we can use a formula established by P. Hagan^. 



Hagan's formula. For a process dF = a{t)c{F)dW , the log-normal implied volatility, at 
leading order in t, is given by 



o-B = a- 



fa 



24 I c \c 



J av j 



0{r) 



with usually fav = {Fq + ^)/2 and a = (l/r /J" a'^{t)dt)^^'^ . 



^Patrick S. Hagan and Diana E. Woodward, Equivalent Black Volatilities, Applied Mathematical Finance 6, 
147-157 (1999). 
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We now consider that for us a can be replaced by o k I ci^j av) ■ Indeed, for small times 



'K 



1 



c(/a.)2 



TcifavY 



^ (\{{AFtf\F,^K) 

lav) JO 

a^{t)c{Ftfdt \Fr^K 



TC-i.fav)'^ Jo 
1 



'K 



Tc{fav)'^ \Jo 
2 



cUavY 



if fav is defined by 



E ( j a\t)c{FtYdt I Fr = Kj = c{favY j a\t)dt. 



Rather than (/o + K)/2, one should take, for small times. 

fav 



K' - fS 



Indeed this is the value that we got for a log-normal model (c(F) ~ F) (see appendix 3). 

Remark. Pierre Henry-Labordere gave a different formula for the implied volatility smile'^ . 

4.3.3 Smiles 

The formulae above give us the exact asymptotic value of the smiles for r = : 



T^O 



= (^T^^(0l,02)) 



1 



(l-/3)(2 + /3) fK-f 
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fa 



with : 

• tan 6i 

• tan 6*2 

• fav = 



yT^/(^./a du/ciu) -p) = ^T^/{u/a ^'^-P^t"^'' ~ p) : 



'°8(t) 



For different values of 5, we get for instance the following smiles (Fig 3.4). 

With all the previous results, one also can compute a first-order (in t) approximation of the smiles, 
which gives for example the following evolution (Fig 3.5). 



•^Pierre Henry-Labordere, A General Asymptotic Implied Volatility for Stochastic Volatility Models, preprint. 



26 CHAPTER 4. APPLICATIONS OF THE PREVIOUS ASYMPTOTICS 




Figure 4.3: Smiles from 6 = 0.5 (down) to S = 0.7 (up). 




Appendix 1 : conditions for non-explosion of 

the volatility 



Let us first remind Feller's Test for Explosions. We consider a process 

dXt = a{Xt)dt + biXt)dWt 

on a segment / of E. We suppose that the functions a and b are Borel-measurable and follow 

the conditions of non-degeneracy and local intcgrability : 



• (Vxe/) (52 (x) >0) 



. (Vx e /) {3e > 0) '-^dy < +^) ; 

The scale function is defined, for a constant c G /, as 

e -I" WTd^. 

One can then define the function 

v{x) = I p'{y) I -TTT^ir^'^y- 
Jc Jc P\z)¥{z) 

The exit time from / is 

S = inf{t ^ I ^ /}. 
We now can write the theorems we need^. 

Proposition 1. If a process Xt follows the conditions above, with Xq ~ x Cz I , then 

• if p{l^) = -oo and p{r^) = +00 then P[S = 00] = P[supo<t<+oo Xt = r] = P[info<t<+oo Xt = 
Z] = 1 ; the process is recurrent; 

• if p{l^) > —00 andp{r^) = +00 then P[supQ^t<+oc Xt < r] = P[\iTat^s Xt = I] = I ; 

• if p{l^) = —00 andp{r^) < +00 then P[supQ^^^^^ Xt > ^] = P[lmit-^s Xt = r] = I ; 

• ifp{l+) > -00 andpir-) < +00 then P[limt^s Xt ^ I] ^ I - P[limt^s Xt = r] = g;_7rff) ; 

Proposition 2 : Feller's Test for Explosions. // a process Xt follows the conditions above, with 
Xq G /, then 

P[S = +00] = 1 v{l+) = v{r-) = +00. 

-■^For a demonstration, the reader can read lOANNls Karatzas, Steven E. Shreve, Brownian Motion and Stochastic 
Calculus, Springer. 
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The case that interests us is 

a{x) = A(/i — x) 



b{x) = u x^ 

with V > Q and 5,X,pi ^ 0. We will apply Feller's test successively for 5^1 and S ^ 1/2, then (5=1 
and (5 = 1/2. 

5.1 If (5 / 1 and (5 / 1/2 

We have 

p{x) = e '^'i'' d^, 



_„ fx A(fi-;)dC 

p'{x) = e 
We then have 

/■^ -2A/( 1-2.5 , A 2-25 /"^ OH r 2A(t 1-2.5 A 2-25 



/(y) 

We will look for equivalents of p or f{y) in the following cases, to apply the proposition 1 or 2. 
If 6 < 1/2 

In this case, one can write 

2dz 2Ati 1-25 A „2-25 

because ^c^'-'ti-^*) '--^(i-*) which is intcgrablc in 0. So we have u(0+) < +oo : by 

proposition 2 P[S < +oo] > 0. 

If (5 > 1 

An analogous calculation as above gives v{0~^) ~ v{+oo) = +oo, so there is no explosion in a finite 
time. However, if A = (no mean reversion) we have p{+oo) < +oo. If A 7^ 0, p(O^) = —00 and 
p{+oo) = +00, so by proposition 1 the process is recurrent. 

If 1/2 < (5 < 1 
We now have 

-2An ^1-25 I A ,,2-25 

/(y) ^ e-"(i-2«) -"(1-*)'^ X csi ^ csi 

y — ^+00 

for y large enough. So w(+oo) = +00. In 0"'", we have, if A > 0, 

-2A£_ 1-25 

So, in this case w(0+) = +00 : P[S < +00] = 0. If A = 0, then 



5.2. IF 5=1 OR 6=1/2 
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which is integrable in 0+. So, in this case, P[S < +oo] > 0. 

Moreover, one can check that for A > 0, we are in the first case of proposition 1, so the process is 
recurrent. 

5.2 If^=lor^=l/2 

Simple calculations (where logarithms appear) show that, for these limit cases : 

• for 5=1 there is no explosion and the process is recurrent ; 

• for (5 = 1/2 there is no explosion if and only if 2\^/v^ > 1 ; in such a case the process is recurrent. 

5.3 Summary 

The cases for which there is no explosion and/or limit for the process are : 

• 5 > 1/2 with mean reversion ; 

• 5 = 1/2 with 2\n/v'^ > 1. 

From a practical point of view, we consider in our note the case 5 € [1/2, 1] with no condition on 
A, /U and v, because our asymptotics are given for short times : the cases of explosions in the interval 
[1/2, 1] would be very pathological. 
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Appendix 2 : 5-geometry 



In this part, we recall some useful results of Riemannian geometry, and we apply it to the special case 
or the upper half plane (y > 0) with the metric 

with S G R"*". We will be particularly interested in the cases S ~ 1 and S = 1/2. 

6.1 Generalities about Riemannian spaces 

In this section, wc note {gij)i<ij<n the coefRcients of a metric G, and we wiU use Einstein's summation 
conventions. 

6.1.1 Christoffel symbols 

The ChristofTel symbols are defined as 

6.1.2 Curvatures 

The Riemann tensor, the Ricci tensor and the scalar curvature are successively defined by 

R = ^9^'^RtJ.,K 

where « ,7 »denotes the partial derivative face to 7. The Gauss curvature is defined in dimension 2 by 

(Rs^tu) ■ V — — i? vol(s, t)vol(w, v), 

that is to say, for ei and 62 orthogonal vectors with norm 1, 

R = (^61,6262) • ei = i?i2^^- 

It coincides with the scalar curvature, in dimension 2. 
In dimension 2, one can show that 



V9 



I ( ^tl) - 1 ( ^n, 



y \ 3 XX / Ox \ Qxx 



with g the determinant of g. 
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6.1.3 Isometries 

Let $ be a function from G to H, with respective metrics G and H. Then $ is said to be an isometry 
iff for every a and b in G 

d{a,b)G = d($(a),$(6))H. 
This global definition is equivalent to the following local condition : for every point Z G G, 

Gz = (V$)5i?$(z)V$z. 

6.1.4 Geodesies 
Minimization of length 

The length of a curve c{t) in a Riemannian space G, between a and b, is defined by the integral 

d{a, b) = J \l < c{t),c{t) >Gdt, 

which is clearly not dependant of the parametrization. 

For sufficiently regular spaces (for instance spaces with scalar curvature always negative, as in 
(5-geometry) , there exists an unique curve j{t) minimizing the length between a and b along 7. This 
length, l{a,b), is the geodesic distance between a and b. A non-trivial property is that this length 
coincides with the minimal energy between a and 6, that is to say the minimum along all ways of 

J <c{t),c{t) >G dt, 

under the normalization condition |c(t)|G = 1- 
The parametric equation 

Let c(t) = [x(r),y(r)] be the position of a point at time r. Then an application of the Euler's 
perturbation formula gives 

V,c = 0. 

In a coordinates system (xi)i<i<„, this can be written, for all k e |1, n] 

ic'' + r'^ix^x^ = 0. 

6.1.5 Killing fields 

A field that preserves the metric is said to be a Killing field. In terms of the Lie derivative, this means 
that 

Cxg = 0, 

so for all fields Y and Z 

.g(VyX,Z)+.g(r,VzX) =0. 
In particular, if the metric is invariant along the direction Xi, then 9^. is a killing field. The previous 
equation then implies, for a way ^(t), g{\/ ■ .dxi,lit)) = 0. If 7 is a geodesic, we then have 

dt9{dx.,j{t)) = .9(V^(,)a,,,7W)+5(a.,,V^^^^7W) 

= + 
dt9idx,nit)) = 0. 
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So, for a given geodesic, one can associate to any Killing field an invariant quantity along the 
geodesic : 

9{dx,,l{t)) = est. 



6.1.6 Jacobi Fields 

A field ^ {t) along a geodesic ^{t) is said to be a Jacobi Field if and only if 

V2^(0+i?(^(t),7)7 = 0, 

where R is the Riemann tensor. This abstract definition becomes more intuitive thanks to the following 
proposition (see Lang^ for a demonstration). 

Proposition. Let 7 : [a, fe]x, [io,^i] ^ G 
be a variation of a geodesic "f{t) such as for every 
s € [a,b] 7(5, .) is a geodesic. Then 



/{t) = dsi{s,t)\ 



s=0 



is a Jacobi Field. Furthermore, to every Jacobi 
Field along 7 one can associate such a variation. 




6.1.7 The Laplace Beltrami operator 

The traditional Laplacian is defined as A(/) = d^^f. In the cuclidean case, this allows to write the 
relation J (pAfdveuci = — J V(y9V/d?Jeuc; for any sufficiently regular function (f> with compact support. 

For any Ricmannian space M, with dimension n and metric M, we want to find an expression of 
the Laplacian that verifies Stoke's theorem, that is to say 



(^A^/dz;^ 



< Wlp, V/ >M dv'' 



where dv^ is the volume element of the Riemannian space (we know^ that it is given by dv^ = 
VdetM dx^...da;".) Let m be the determinant of the metric M. As we have (we use Einstein's 
summation convention) 

J < V(/3, V/ >M dwM = J < V(/3, V/ >M da;-' . . .dec" 

9''^§^V^dx\..dx- 

Ox3 OX^ 

j < Vv?, V/ >M d.T^ . . .dx" = -y"v3--^^- w*^^^ V^da;^...dx", 

it is necessary and sufficient to define the Laplace Beltrami operator as 



"'^Serge Lang, Fundamentals of Differential Geometry, Springer. 

^JURGEN JOST, Riemannian Geometry and Geometric Analysis, Springer Univcrsitext, 2002. 
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One of the advantages of this intrinsic formulation of the 
Laplacian is that it is invariant by isometry. More precisely, 
if $ is an isometry between two spaces G and H, then for 
every function / from H to R this diagram is commutative : 

A^(/o$) = A^(/)o$. 



\ I / 



Indeed, as an isometry conserves the scalar product (so the 

volume element, too), we have for every sufficiently regular function cp from G to R with compact 
support. 



A^(/o$) if d/ 



J <V(/o$),V^ >G dv^ 
J < V/,V(v?o$-i) >H dv^- 
J A«(/) ^o$-i dv^ 
J AH(/) o $ ^ d«« 



As if is any sufficiently regular function with compact support, this gives the commutation relationship. 



6.2 A special case : the (5-geometry 

The 5-geometry is defined on the upper half plane {y > 0) with the metric 



2 da;^ + dy"^ 



ds' = 



with (5 e R+. We will represent it with the matrices 



,,2(5 



G 







v"' 





6.2.1 The geodesies 

After some calculation, the only Christoffel symbols that are not equal to are 

. Vly^\l2g^y{dyayy) = -^ly- 
. = l/25^n-5y5x.) = Sly. 

Parametric equations 

The parametric equation for geodesies is here 

/ i = 

1 y = |(y -•^') ■ 
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Another (equivalent but more practical) differential system is the following, 

2 . ■ 2 

= ctel 

y 

= cte2 

This system needs some explanation : 

• the first equation expresses that 5(7,7) is constant ; 

• the second one expresses that d/dx is a killing field : the metric coefficients do not depend on x. 

With these equations we can get the graph of these geodesies, leaving from [0, 1] with a speed [1, 0] 
and for r e [0, 100] (Fig 5.1 and 5.2). 




0.2 0.4 0.6 0.8 1 



Figure 6.1: ^ — 1 




1 2 3 4 5 



Figure 6.2: 6 = 1/2 



Explicit equation 

Eliminating dt in the second system of the previous paragraph, we get (imposing a maximal possible 
y equal to 1) 

y^dy 



dx = 

v/1 - y^' 

An integration gives the following expression for x{y), imposing a;(0) = : 
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get the following graphic for different values of S (Fig 5.3). 




0.2 0.4 0.60.8 1 

Figure 6.3: 6 G [0.3, 2] 

An interesting point is that for 6 — 1/n (n € N), Mathematica always finds an explicit formula. 



Explicit formulae for x{y) 


S 




1 




1/2 


arcsin ^/y- \/y- y"^ 


1/3 




1/4 


\ (3 arcsin _ _ ^(3 + 2^)y^/^) 
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6.2.2 Curvatures 



For our metric, the only non-zero curvatures are 



yxy 

TiX 

yy^ 

Ry 



Ry 



xxy 
xyx 



^5 



<5 



''"yy 



R = y 



2sd 5 _ 



dyy 



6.2.3 Homogeneity 

Notice first that if 7{t)tg<t<ti ^ ^^^'^ '^^ ^^'^ upper half plane (not necessarily a geodesic), then one 
can easily know the length of A7 (A £ M+) : 



/(A7) 



g{X-f,Xj)dt 



.2 .2 

Xx + Xy 



dt 



= X 



25 



dt 



Z(A7) = X'-'l{j) 



Another important homogeneity property is that if 7 is a geodesic, so is A7 : replacing x and y in 
the parametric equations defining the geodesies by Ax and Xy keeps the equalities true. 



6.2.4 Distance from a point to a straight line 



For our metric, one can easily check that vertical lines are geodesies, so the distance from a point to 
a vertical is gotten when the geodesic is perpendicular to the vertical axis. 

So the homogeneity properties make it possible to calculate this length. 
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A 




We have a point M = {X,Y). If the smallest geodesic is a « reference » geodesic then one can 
calculate the abscisse x of the intersection of on the schema. Then 

6.3 The first-order expansion at the origin. 

To give Ki at the origin {xQ,yo), we need to evaluate the terms in (2.4). In this equation g can be 
replaced by \J djZid) = e^h. First notice that {,^\xo,yQ) — 1- Moreover, from the differential 
equation in the corollary 1, we have Z(d) = c? + — ^^^|p^c?'^ + 0(d^) so 



li?(xo,2/o)x^+2/^ .0(d3) (2.1) 



YZ(d) y d+^iiMd3 + o(rf4) 2 3! yf 

for near points (a;o,yo) and (xo + a;,?/o + v)- We also need an expansion of W up to the third order. 
This is quite technical but necessary. We consider the points a = (xq, yo) and h = (xq + x, j/q + ?/)j and 
we note VF(a, 6) = iy(x, y). Then a classical Taylor expansion gives 

W{x, 2/) = VKo + (9,iy)o X + (dyW)^ (d^^W)^ + 1 (a,,W^)o y^ + {d^yW)^ xy + 0{d^). 

The evaluation of the partial derivatives with respect to y is straightforward, because we know that 
the vertical hnes are geodesies for the 5— space (cf appendix 2). More precisely. 



WiO,y)^y 



' fy{0,ty)dt 
ivo+ty)^^' 



and 

1 fr^ My)o 



yo 

The evaluation of the partial derivatives with respect to x is more complicated, because the geodesic 
between (xq, yo) and (xq + x, j/o) is not a straight line anymore. Anyway, we will show that the calcu- 
lation is the same as if it was a straight line indeed ; to do so, we will use the symmetry with respect 
to the line of points with abscise xq + x/2 (cf appendix 2). 
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A simple look at the geodesies in 5-geometry shows that the following schema is true : if a: is the 
length of the straight line, then the length of the geodesic is 0{x) and the height is 0{x^). 




(xo + X, yo) 



Let us Wi (resp W2) be the work of the field / along the geodesic (resp the straight line). We now 
show that 

Wi =W2 + Oix^). 

The key point is that the geodesic has a symmetry axis, which implies tan0d/2 = 0. This allows 
us to write 

W1-W2 = 



(/i-/2)-dZi+ J f2-idh~dh) 

0(a;=) 

0{x')+ I [(/Do + (9./|)o^2 + 0(/2)] tanM/2 



0{xn+ I {d^fDo^tanedh 

0{x) 0{x) 



W1-W2 = Oix^) 



Thus we have 



idxxWi)o = id^xW2)o, 

which can be easily calculated, like for the partial derivatives with respect to y. We finally get 

{dxxWi)o = -^{d^fx)o- 
Vo 

We now can evaluate (2.4). We have h{xQ,yo) ~ 1, so the equation (2.4) at the point {xo,yo) 
becomes 

(/u)o=[(V + /).g]^, 

that is to say, for g — h, 
1 



yf [e^A + 2Ve^ • V/i + hAe'^h] „ + [e"^ fh + hfe"^] ^ 



25 



(A/i)o +2(Ve^)o • (V/i)o + \VW\l + (AM^)o 

(1) (2) (3) (4) . 



+ {fh)o + {fW)o 

(5) (6) 
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We now have to evaluate each term : 

• from (2.5) one can easily conclude that (1) = R{xQ,yo)/ {Si/q^) and (2) = (5) = ; 

• the first-order approxhnation for W gives (3) = |/o| /uo^^ and (6) = |/o| /vo^^ ; 

• the second-order evaluation of W gives 



(4) 



1 



yo 



25 



(div/)o 



yo 



We then finally have 



(^i)o 



25-2 



(div/)o {fy)o , 3|/op 



yo 



2 yr 



6.4 Two new functions 

Thanks to the above homogeneity properties, we only need to calculate the distances and the Jacobi 
fields on the standard geodesic (symmetric with respect to the ordinatcs axis and y[0] ~ 1). Here is 
the natural way to do all the necessary calculations on the standard geodesic : 

• the function Sind gives the position that minimizes the distance from a point to a line ; 

• the distance function is specially easy to compute on the standard geodesic : we have a formula 
with only the coordinates of the end points ; 

• then, to compute the Jacobi fields,we use the fact that it is a 2-dimensional linear space. 



6.4.1 The function sind 

We write a;(y)the function such as 




Geometric definition of function sind 

The intuitive definition of function sind is given by the following drawing. 

y 




(e,0) 

Figure 6.4: The geometric view of function sind 
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Analytical definition of function sind 

For X > 0, we have 



dx —y 



2S 



so 

1 

yi = (sin^i)^ . 

A simple calculation then gives 



if f < 01 < TT, -x(yl) - ^ 



2 — -'i ^ " 7 -^KU^J tan Si 

Then the function sind((5, 9i, 62) is the unique solution of the following equation 
f if e + -V > 0, e + = x(smd) 

} tan $2 — ' tan ^2 ^ 



if e + -V < 0, e + = -xfsind) 



(2.2) 



Diff"erential properties of function sind 



The equation (1) defines sind in an implicit way. One can now get the differentials of sind, differentiating 
(1). The equations should be different in 4 zones of [0,7r] x [0,7r], but we just have 2 such zones in 
reality. 

if e + -V > 0. ( , + -Vl ae.sind = (i - l) (sin^^^-^ 

tan 92 - VVl-Sind'-* '^""V ;V IJ 

ife+tl^>0, f ^MM;= + a,,sind = 



Vi-sind=* tane. ;-''2 - (,i„e,) 



ife+-V<0, ( - , ^i"d^ + -V ) c>g,sind ^ 

tane2 — ' I Vl-Sind tan 62 / ''^ (sin6l2) 

Examples of function sind 

One can check that for (5=1 (geodesies are circles) the sind(0i,6'2 — sin(0i)/ sin(02))- 

The undefined values appearing on the frontier above is due to the difficulties Mathematica has to 
solve the equation (2.1) for 

e+ « 0. 

tan 02 

6.4.2 The distance 

Thanks to the homogeneity properties, we only need to know the distance between the point (0,1) and 
any other point of the standard geodesic. This distance is exactly 



d((0,l)^ (xi,2;i)) = 



' da;^ + dy^ 
^23 



dy 



d((0,l)^ (xi,yi)) 



x,-yl-'./T^ 
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Figure 6.5: S — 1 




Figure 6.6: 6 = 1/2 



thanks to links between hypergeometric functions. 
6.4.3 The Jacobi field 

The equation giving the Jacobi field on any geodesic is linear, so it is sufficient to evaluate it on the 
standard geodesic. In order to evaluate the Jacobi field on the standard geodesic with any initial point 
and initial increase, we just need to calculate two specific geodesies and then use linear combinations. 
The most natural initial conditions are : 

• Z{0,1) = and Z{0,1) = 1 ; 

• Z{0,1) = 1 and Z{0,1) = 0; 

The differential equation giving the evolution of a Jacobi field (we just need to know its length 
because it remains orthogonal to the geodesic) is 

Z + K Z = 
^ Z - 6y^^-^Z = 0. 
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For 5 = 1, the equation above has the evident solution Z{d) = Acoshd + Bs'mhd. For S ^ 1, we 
need to reformulate the equation above in terms of the coordinate y. We get 



(1 - y'^)dyyZ + S —dyZ -—Z = 0. 



The problem is that this equation is not lipschitzian for y = I. This is the reason why, to get a 
great approximation of the solution, we use the initial conditions 

Zi[l-e] = r Z2[l-e] = 1 
Zi[l-e] = 1 ' I Z2[l-e] = 

with e small enough (typically e = 1/100). If we take the x-coordinates, we then just need to take : 

• if X > 0, Zi[x] := Zi[y[x]], if x < 0, Zi[x] := -Zi[y[x]] ; 

• if a; > 0, Z2[x] Z2[y[x]], if x < 0, := Z2[y[x]]. 
Then, linear combinations of Zi and Z2 give any Jacobi Field. 
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6.5 Variation formulas for the (5-model 

The reader who does not know anything about the foUowing variation calculations may read the famous 
books by Milnor^ or Lang"*. 



6.5.1 The first variation formula : minimization implies orthogonality 

We consider sufficient conditions of regularity 
(uniqueness of the geodesic between two points 
etc) . Then the distance from a point Zq to a line 
is minimal in a point Zmm (with vector along the 
geodesic 'y{Zmin) and along the line directed by 
v{Zrmn)) if and only if 



< liZ,mn),v{Z,nm) >H= 0. 




As H is locally euclidcan, this is equivalent to 
the usual orthogonality. Let us justify this rela- 
tionship. For a variation 7(5, t) (where t repre- 
sents the proper time, that is to say the length), 

if 7(s, •) is a geodesic for every s, then the length Figure 6.7: Geodesies between a point and a hy- 
and the energy are equal, that is to say persurface in a curved space : orthogonality at the 

critical point. 



d{s) = V< dais, t), dais, t) >Mdt = 2 y < dais, t), dais, t) >H At. 

This is a simple consequence of the definition of the proper time : < dais, t), dais, t) >h= 1- We then 
can write 



dsdis) 



1 



ds < dais, t),dais, t) >M dt 



-y (< Va.dais,t),dais,t) >n + < dais,t),Va^dais,t) >M)dt 
\ j i< ^d,dais,t),dais,t) >H + < dais,t),\Id,d.ais,t) >u)dt 
< "^dfdslis, t),dais, t) >H dt 

idt < dsjis,t),dais,t) >M~ < ds-fis,t),\/otdais,t) >M)dt 
dt < dsjis,t),dais,t) >H dt 

= < dais, t), dais, t)>M\lzo 
dsdis) = < dais, t), dais, t) >M r=''. 

The calculation above requires some explanation : 

• (2.1) : we use the well-known relation dx < u,v >=< Va^it, v > + < u, Va^w > ; 

• (2.2) : the connection V is torsion-free ; 



(2.3) 
(2.4) 

(2.5) 
(2.6) 

(2.7) 



^J. MiLNOR, Morse Theory, Annals of mathematics studies, Princeton university press. 
''Serge Lang, Fundamentals of Differential Geometry, Springer. 
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• (2.3) 

• (2.4) 

• (2.5) 



we also use the relationship dx < u,v >=< Vg^u, w > + < u, Vg^w > ; 
7 is a geodesic for every s, so \7 g^dt'^{s,t) = 0; 

the initial point docs not move with the variation, so ds"f{s, i)|t=o — 0. 



6.5.2 Second variation for any strike 

In the following, we will have to calculate partial derivatives with respect to s, so for simplicity we will 
write rj = dg. We will then often use the following expression for the connection with direction 77, 



dgUx 




dsUy 


+ 



KyVxUy + Tl^r^yUx 

^xxVxUx + ^yyVy'^y 





dgUx 


5 \ 
+ - 




dgUy 


V 



— cosoiUy — smtllUx 
cosOiUx — sinOiUy 



(2.8) 



where 9i is the angle between the straight line and the abscises axis. 

To calculate the second derivative, we derive it from the expression (2.5), 

d,sd{s) = <d,-f{s,t),dtjis,t)>M\'=^ 

= < V9Aiis,t),dais,t) >H 1*="^+ < a,7(s,t), Va, 9*7(5, i) >H 1*='' 
= < Vs,<9,7(s,i), 9*7(5, t) >H r='^+ < 9.7(s,i), Va,9,7(s,f) >h 1*=^^ 



1 



< V„r/, dais, t) >H 1*-'' + -dt < V, 1] >H I 



\t=d 



(2.9) 



(1) 



(2) 



where we use that the connection is torsion-frcc, and the relationship of derivation of the scalar product. 
Let us evaluate each of both terms. 

• We first have to calculate the covariant derivative 



5 

y 



-2 cos 6*1 sin6'i 
"2 fl- sin^ 61 



cos (71 

at the point where 9i is the angle between the abscises axis and D. Moreover, 

7 = ±/ 



sm Ui 
- cos 61 



where the sign + or - depends on which side the point Zq is with respect to D. So we have 

(1) = ±775+1 

mm 

• As parallel transport keeps orthogonality, we have ?7|*=° = ■§^^, where ~u is the orthogonal of 
the speed vector with norm 1 (no matter its direction). So we have 

Zid) 



(2) 



To sum up, we have 



dy' 



rd(Y,r. 



Z{d)Y?A 
( Z{d) 



h \Z{d)Y^ 



2r± 775+1 cos 01 

min ^ min , 



where the sign depends on the following configurations. 
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6.5.3 Second, third and fourth variations at the money 

We will also need the second, third and fourth derivatives, for d — > 0. 

• we have $" = dd" ^ l .^^ ; 

• $(3) = ^^(3) 

SO we have to calculate the third derivative of the distance. Let us derive it from 
the formula (2.7) : 

dsssd{s) = < \7n'^rjr],dt'-f{s,t) >h +2 < V,,7?, V^7(vi) >h + < yr,V,)<9t7(s, t) , V >h ■ 

(1) (2) (3) 

We have dx (1) and dx (2) n. The term (3) is more difficult to evaluate : it represents 

the variation with s of the Jacobi-field increase, which is clearly maximum for s = ; this is the 
reason why the only term that remains in the calculation of (3) is the Christoffel part of the 
covariant derivative. A simple calculation then gives 

$(3)|/^=/o -3(5 



Yo^'+' sin2 0, ■ 



Unfortunately, at this time, we don't have properly derived a good expression for However, 
an analogy with the 6=1 case makes us believe that 

Anyway, this expression is exact for 6=1. 



Appendix 3 : calculation of fav for" the 
implied volatility formula 



We describe hereafter the derivation of fav introduced in [Hagan [10]] as wc need it in chapter 4.3. 

We consider the process dF — a{t)FdW {F{0) = /q), and wc want to calculate the value fa 
defined by 



E 



^ a''{t)Ft^dt \Fr^K^^ fav^ ^ a^{t)dt. 



7.1 The solution of the process 

We know that Ft is a martingale, with explicit value 

i?, = /oe^o"(^)div.-i/>^(.)d._ 

This is the reason why 



£ a'F^dt \Fr = K^ = E a' f^e^^ -du^-Jo' c-ds^^ | ^^^j; adw-i Jo^ "^d. ^ ^ 

= fie- «'d.E (^^2 /„' adw^^ | £ ^ 



E 



with c = ln(A'//o) + 1/2 a^ds. We now just need to know the conditional law of adiy with 
respect to adW = c. This is a problem that evokes the brownian bridge, as seen in next part. 

7.2 An equivalent to the brownian bridge 

Lemma The law of {J^ adW \ adW — c) is the same as the law of 

adW- -[Q / adW+ I" „^ c. 

Jg a'^ds Jq Jg a''ds 

Proof. We follow the traditional proof of the brownian bridge. 
• first of all, J^^ adW and adW — |° "^^^ JJ" adW are independant, because : 

— {J" Q;diy)„>o is a gaussian process, so (/„ adW, L adW — " C adW) is a gaussian 

Jo ^ 

vector ; 
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— so its coordinates are independant iff their covariance is zero ; this is true (the coefficient 
{° 1^ has been chosen for it) ; 

• let us (phc a sufficiently regular and bounded function. Wc know that, if U and V are independant 
processes, then E {ip{U, V) \ V = v) = E {ip[U, v)). So we can write 

E (^j^ adW^ I ^ adW = c 

/ / /•* r'a^ds r f*aMs r \ T \ 

= E / adW - — / adW + ■'° , / adW \ / adW = c 

\ \Jo Jo " ds Jo Jo a ds Jo J Jo J 

This means exactly the result of the lemma : both processes have the same law. □ 



7.3 Final calculation 



We now can directly calculate the expectation (here /3 is the suitable function depending on s < t or 

s > t) 



E ( e^/o^dH'dt I j adW^c 



E e 



r* Q^d 



2 ( / ' adW- lifL^ / - adw+ HfL^ 



dt 



This is how we get the final expression for /, 



ft „2(^g 

''ir:^s\ks;i3\s)ds 



fav — 






dt 






/J" "^ds 



One can check, after tedious calculations, that 



fav 



^-°\log(f 



If one does not want such a complicated expression, it can be approximated, for K 
(/o + K)/2. 



fo, by 



